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Abstract 

We discuss the derivation of the path integral representation over gauge degrees of freedom 
for Wilson loops in SU{N) gauge theory and 4-dimensional Euclidean space-time by using 
well-known properties of group characters. A discretized form of the path integral is naturally 
provided by the properties of group characters and does not need any artificial regularization. 
We show that the path integral over gauge degrees of freedom for Wilson loops derived 
by Diakonov and Petrov (Phys. Lett. B224 (1989) 131) by using a special regularization 
is erroneous and predicts zero for the Wilson loop. This property is obtained by direct 
evaluation of path integrals for Wilson loops defined for pure SU (2) gauge fields and Z{2) 
center vortices with spatial azimuthal symmetry. Further we show that both derivations 
given by Diakonov and Petrov for their regularized path integral, if done correctly, predict 
also zero for Wilson loops. Therefore, the application of their path integral representation of 
Wilson loops cannot give a new way to check confinement in lattice as has been declared by 
Diakonov and Petrov (Phys. Lett. B242 (1990) 425). From the path integral representation 
which we consider we conclude that no new non-Abelian Stokes theorem can exist for Wilson 
loops except the old-fashioned one derived by means of the path-ordering procedure. 
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1 Introduction 



The hypothesis of quark confinement, bridging the hypothesis of the existence of quarks and 
the failure of the detection of quarks as isolated objects, is a challenge for QCD. As a criterion 
of colour confinement in QCD, Wilson [1] suggested to consider the average value of an operator 

W{C) = ^i,Vce'9 §cdx,A,{x) ^ ^tr[/(C7..), (l-l) 

defined on an closed loop C, where A^{x) = is a gauge field, (a = 1, . . . , A^^ — 1) 

are the generators of the SU{N) gauge group in fundamental representation normalized by the 
condition tr {t"-t^) = 5'^''/2, g is the gauge coupling constant and Vc is the operator ordering 
colour matrices along the path C. The trace in Eq.(|l.l|) is computed over colour indices. The 
operator 

makes a parallel transport along the path Cyx from x to y. For Wilson loops the contour C 
defines a closed path Cxx- For determinations of the parallel transport operator U{Cyx) the 
action of the path-ordering operator Vcy^ is defined by the following limiting procedure [2] 

U{Cyx) = VcJ' ^"^^^^"^ = Jim n U{Cx,x,.,) = 

k=l 

n 

= limUiCyx„^,) ...UiCx,x,)U{Cx,x)= lim T\ 9 i^k - Xk-i) ■ A{xk-i) ^ (1,3) 

n — ^oo n — ^oo 

k=l 

where Cx,.x,,_i is an infinitesimal segment of the path Cyx with xq = x and Xn = y- The parallel 
transport operator U{Cx,.xu-i) infinitesimal segment Cx,.xk-i is defined by [2]: 

U{Cx,x,^,) = -^^^^^^-1 dz^A^iz) ^ ^ig^^^ _ ^^_^) . A{xk_^)_ (-^ ^^ 

In accordance with the definition of the path-ordering procedure (|1.3| ) the parallel transport 
operator U{Cyx) has the property 

U{Cyx)=U{Cyx,)U{Cx,x), (l-S) 

where xi belongs to the path Cyx- Under gauge transformations with a gauge function 

A^{z) ^ A^'^iz) = n{z)A^{z)n\z) + i d^^{z)n\z), (1.6) 
the operator U{Cyx) has a very simple transformation law 

U{Cyx) ^ U^'iCyx) = n{y) U{Cyx) ^^x) . (l.7) 

We would like to stress that this equation is valid even if the gauge functions Vt{x) and ^{y) 
differ significantly for adjacent points x and y. 

As has been postulated by Wilson [1] the average value of the Wilson loop < W{C) > in the 
confinement regime should show area-law falloff [1] 

< VF(C) >~ e"^-^, (1.8) 
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where a and A are the string tension and the minimal area of the loop, respectively. As usually 
the minimal area is a rectangle of size L x T. In this case the exponent a A can be represented 
in the equivalent form a A = ViJ-^ T, where V{L) = aL is the interquark potential and L is the 
relative distance between quark and anti-quark. 

The paper is organized as follows. In Sect. 2 we discuss the path integral representation for 
Wilson loops by using well-known properties of group characters. The discretized form of this 
path integral is naturally provided by properties of group characters and does not need any arti- 
ficial regularization. We derive a closed expression for Wilson loops in irreducible representation 
j of SU{2). In Sect. 3 we extend the path integral representation to the gauge group SU{N). 
As an example, we give an explicit representation for Wilson loops in the fundamental repre- 
sentation of SU{3). In Sects. 4 and 5 we evaluate the path integral for Wilson loops, suggested 
in Ref.[3], for two specific gauge field configurations (i) a pure gauge field in the fundamental 
representation of SU{2) and (ii) Z{2) center vortices with spatial azimuthal symmetry, respec- 
tively. We show that this path integral representation fails to describe the original Wilson loop 
for both cases. In Sect. 6 we show that the regularized evolution operator in Ref.[3] representing 
Wilson loops in the form of the path integral over gauge degrees of freedom has been evaluated 
incorrectly by Diakonov and Petrov. The correct value for the evolution operator is zero. This 
result agrees with those obtained in Sects. 4 and 5. In Sect. 7 we criticize the removal of the 
oscillating factor from the evolution operator suggested in Ref . [3] via a shift of energy levels of 
the axial-symmetric top. We show that such a removal is prohibited. It leads to a change of 
symmetry of the starting system from SU{2) to U{2). Keeping the oscillating factor one gets 
a vanishing value of Wilson loops in agreement with our results in Sects. 4, 5 and 6. In the 
Appendix we evaluate the coefficients of the expansion used for the path integrals in Sects. 4 
and 5. 

2 Path integral representation for Wilson loops 

Attempts to derive a path integral representation for Wilson loops ( |1.1| ) , where the path ordering 
operator is replaced by a path integral, have been undertaken in Refs.[3-5]. The path integral 
representations have been derived for Wilson loops in terms of gauge degrees of freedom (bosonic 
variables) [3,4] and fermionic degrees of freedom (Grassmann variables) [5]. For the derivation 
of the quoted path integral representations for Wilson loop different mathematical machineries 
have been used. Below we discuss the derivation of the path integral representation for Wilson 
loops in terms of gauge degrees of freedom by using well-known properties of group characters. 
In this case a discretized form of path integrals is naturally provided by the properties of group 
characters and the completeness condition of gauge functions. It coincides with the standard 
discretization of Feynman path integrals [6] and does not need any artificial regularization. 

We argue that the path integral representation for Wilson loops suggested by Diakonov and 
Petrov in Ref. [3] is erroneous. For the derivation of this path integral representation Diakonov 
and Petrov have used a special regularization drawing an analogy with an axial-symmetric top. 
The moments of inertia of this top are taken finally to zero. As we show below this path integral 
amounts to zero for Wilson loops defined for SU{2). Therefore, it is not a surprise that the 
application of this erroneous path integral representation to the evaluation of the average value 
of Wilson loops has led to the conclusion that for large loops the area-law falloff is present for 
colour charges taken in any irreducible representation r of SU{N) [7]. This statement has not 
been supported by numerical simulations within lattice QCD [8]. As has been verified, e.g. in 
Ref. [8] for SU{3), in the confined phase and at large distances, colour charges with non-zero 
A^-ality have string tensions of the corresponding fundamental representation, whereas colour 
charges with zero A^-ality are screened by gluons and cannot form a string at large distances. 
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Hence, the results obtained in Ref.[6] cannot give a new way to check confinement in lattice as 
has been declared by Diakonov and Petrov. 

For the derivation of Wilson loops in the form of a path integral over gauge degrees of 
freedom by using well-known properties of group characters it is convenient to represent W{C) 
in terms of characters of irreducible representations of SU{N) [9-11] 

Wr{C) = ^^X[Ur{C..)], (2.1) 

where the matrix Ur{Cxx) realizes an irreducible and d^-dimensional matrix representation r of 
the group SU{N) with the character x\Ur{Cxx)] = ^"^WriCxx)]- 

In order to introduce the path integral over gauge degrees of freedom we suggest to use 

j D^rXlUM] xi^rVr] = J xWrVr], (2.2) 

where the matrices Ur and Vr belong to the irreducible representation r, and D^^- is the Haar 
measure normalized to unity J DO,,. = 1. The completeness condition for gauge functions fi^ 
reads 

By using the completeness condition it is convenient to represent the Wilson loop in the form 
of the integral 

Wr{C) = ^ [ Dnr{x)x[nr{x)Ur{CxxH{x)]. (2.4) 



dr 

According to Eq. (|1.3| ) and Eq.( |1.5D the matrix Ur{Cxx) can be decomposed in 

Ur{Cxx) = lim Ur{Cxxn-i)^r{Cxn-ix„-2) ■ ■ ■^r{Cx2Xi)Ur{Cxix)- (2.5) 



Substituting Eq.( |2.5D in Eq.(^jJ) and applying (n — l)-times Eq.(2^) we end up with 

Wr{C) = ^ Y\m^ j . . . j DVLr{xi) . . .Q.r{Xn)drX[^^ 

. . . drX[^r{xi)Ur{Cx,xM{Xn)]- (2.6) 

Using relations 9.r{xk)Ur{Cxf,xu-i)^l{xk-i) = U^iCx^xk-i) we get 

WriC) = 1 • • • / DQrixi) . . . DnriXn) drXp^ {Cx„x„_,)] ■ ■ ■ drX[U^ {Cx.xJ]- (2.7) 

The integrations over Qr{xk) (k = 1, . . . ,n) are well defined. These are standard integrations 
on the compact Lie group SU{N). 

We should emphasize that the integrations over J7r(xfc) (/c = 1, . . . , n) are not correlated and 
should be carried out independently. 

Since Eq.( |2.3| ) is the completeness condition for group elements, the discretization of Wilson 
loops given by Eqs. (|2.6|) and ( |2.7D reproduces the standard discretization of Feynman path 
integrals [6] where infinitesimal time steps can be described by a classical motion. Therefore, 
the discretized expression (2.7) can be represented formally by 

^riC) = 4 / n [drDQrix)] Xp^iCxx)]- (2.8) 
x&C 
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Conversely the evaluation of this path integral corresponds to the discretization given by Eqs.(|2^ 
and ( |2.7| ). The measure of the integration over i^r(x) is well defined and normalized to unity 

J Y[ DQrix) = Imi^ J DQr{Xn) j DVLr{Xn-l) ■■■ j D9.r{xi) = 1. (2.9) 

Thus, for the determination of the path integral over gauge degrees of freedom (|2.8|) we do not 



need to use any regularization, since the discretization given by Eqs.(2^) and (^J) are well 
defined. 

We would like to emphasize that Eg. ( ^.8]) is a continuum analogy of the lattice version of 
the path integral over gauge degrees of freedom for Wilson loops used in Eq.(2.13) of Ref.[ll] 
for the evaluation of the average value of Wilson loops in connection with Z{2) center vortices. 

Now let us to proceed to the evaluation of the characters xW^iCx^Xk-i)]- Due to the 
infinitesimality of the segments Cx,.x,^_i we can omit the path ordering operator in the definition 
of U^{CxkXk-i) [2]- This allows us to evaluate the character x[U^{CxkX^_J\ with U^{Cx^Xk-i) 
taken in the form [2] 

U^{Cx,x,.,) = expig f dXf,A)i{x). (2.10) 



k^k-1 



Of course, the relation given by Eq.( 2.10| ) is only defined in the sense of a meanvalue over an 
infinitesimal segment Cxf.Xk-i- Therefore, it can be regarded to some extent as a smoothness 
condition. Unlike the smoothness condition used by Diakonov and Petrov [3] Eg. ( ^.10 ) does not 
corrupt the Wilson loop represented by the path integral over the gauge degrees of freedom. 

The evaluation of the characters of U^{Cxf,Xk-i) given by Eq.( 2.1C| ) runs as follows. First let 
us consider the simplest case, the SU(2) gauge group, where we have r = j = 0, 1/2, 1, . . . and 
dj = 2j + 1. The character xiU-^iCx^Xk-i)] is equal to [9,10,12] 

x[U^{Cx,x,^,)] = E <m,\U^{Cx,x,^,)\m,>= 

mj=-j 

= J2 e^"^J'^[^^fc^fe-i'^^], (2.11) 

mj=-j 

where ruj is the magnetic colour quantum number, \mj > and rrij ^[Cx^.Xk-i': the eigen- 

states and eigenvalues of the operator 

^Cx,x,.,;A''] = g dx^A^ix), (2.12) 



i.e. ^[Cxf,xk-i\ A^W^j >= ''^j^[Cxf,Xh-i\A^]\^j{xk-i) >■ The standard procedure for the 
evaluation of the eigenvalues gives ^[Cxf.Xk-iJ the form 

^[Cx,x,.,;A^] = g j ^g^,[A^]{x)dx^dx,, (2.13) 

where the metric tensor can be given formally by the expression 

g,M'']{x) = '^^^KA^]{x). (2.14) 

In order to find an explicit expression for the metric tensor we should fix a gauge. As an example 
let us take the Fock-Schwinger gauge 

x^A^{x) = 0. (2.15) 
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In this case the gauge field A^{x) can be expressed m terms of the field strength tensor G^iy{x) 
as follows 

1 

A^{x) = J dssXaGaf_i{xs). (2.16) 



This can be proven by using the obvious relation 

XaGa^{x) = XadaA^[x) — Xad^A^^x) — ig[XaAQ,[x) , A^(^x)] = 

d_ 

' dXa 

valid for the Fock-Schwinger gauge XaAa{x) = 0. Replacing x ^ xs we can represent the r.h.s. 
of Eq.( p.l7 ) as a total derivative with respect to s 

d ^ , _ d 
' dxn ^ ds ' 



-Af,{x), (2.17) 



sXa_Gc^,{xs) = Af,{xs) + Xa-^ — A^{xs) = —[sA^{xs)]. (2.18) 



Integrating out s E [0, 1] we arrive at Eq. (|2.16| ). 

Using Eq.( p.l6| ) we obtain the metric tensor gf^i,[A^]{x) in the form 

1 1 

g^^[A^]{x) = 2xaXi3 J J dsds'ss'ti[G^^{xs)G^^{xs' )] = 





1 1 



2xc,x^ j J dsds'ss'tr[n{xs)Gaf,{xs)n\xs)n{xs')G(3^{xs')n''{xs')]. (2.19) 







For the derivation of Eq.( p.l9|) we define the operator ^[Cx^.Xi._i', A^] of Eq.( 2.12| ) following the 



definition of the phase of the parallel transport operator U{Gxf.Xk-i) given by Eq.( |1.4| ) [2] 

^Gx^xk_,; A^] = g[ dx^A^ix) = {xk- Xk-i)f,A)lixk-i) = 



1 



= {xk-Xk-i)^, j dssxt_^^G^^{xk-is). (2.20) 



The parameter s is to some extent an order parameter distinguishing the gauge functions ^{x^) 
and 0(xfc_i) enterin g the relation 9.{xk)U {Gx^Xk-i)^"^ {xk-i) = U^iCx^Xk-i)- 

Substituting Eq.(|2]T|) in Eq.(^ we arrive at the expression for Wilson loops defined for 
SU{2) 

Wi(G) = 3— :r lim 

^ (2i + l)2n^oo 



/ on, (X.) (2j + 1) ± e'' "^5"^ ^^^1-^ ^9,AA^]{x)dx,dx, 



(n) 

m =—7 



J Dfljixn-i) {2j + 1) e 



^ igmf ^^/c,„,^_j ^Jgf,^[A^]{x)dXf,dx^ 



(n-l) 

m. =—7 



^ i g mf^ Jc^^^^ ^g^^[A^]{x)dx^dxy 



DQj{xi){2j + 1) e^-j JC^,., (2.21 



(1) 

m\ =—7 
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The magnetic quantum number ttij ^ (/c — 1, . . . , 71) belongs to the infinitesimal segment Cx/^^ix^.^ 
where Cx„+,x„ = C^^x^. 

In compact form Eq.( 2.2l| ) can be written as a path integral over gauge functions 



(2i + 1) J {mj{x)} 

The integrals along the infinitesimal segments Cxi^xh_i we determine as [2] 

^ Crr 



mj{x) ^g^y[A^]{x) dx^dx^ = mj{xk-i) gfj,u[A^]{xk-i) Ax^ Ax^ 



mf ^'^Jgf,,[A^]{xk^i)Ax^ Ax,. (2.23) 



where Ax = Xk — Xfc_i. 

Comparing the path integral ( 2.22| ) with that suggested in Eq.(23) of Ref.[3] one finds rather 
strong disagreement. First, this concerns the contribution of different states rrij of the repre- 
sentation j. In the case of the path integral ( |2.22| ) there is a summation over all values of the 
magnetic colour quantum number mj, whereas the representation of Ref.[3] contains only one 
term with rrij = j. Second, Ref.[3] claims that in the integrand of their path integral the exponent 
should depend only on the gauge field projected onto the third axis in colour space. However, 
this is only possible if the gauge functions are slowly varying with x, i.e. Q{xk)^^ {xk-i) — 1. 
In this case the parallel transport operator U^{Cxf,Xk-i) would read [13] 

U^iCx^x^_,)=expig f dx^A^{x) = l + ig{xi-Xi_i)-A^{xi-i), (2.24) 

and the evaluation of the character x[Uf'{Cx^.x^.^l] would run as follows 



< mj\[Uf{Cx,x,.^)]\mj >= 1 + {t'^U^n.Agixk - Xk-i) ■ [A^ixk^i)]^"^ = 

o iqfn dxnmi(x)\A^(x)]^^^ 

= l + m,-i5(^fc-Xfc_i)-[A^(x,_i)](3) = e ^•^'^^^^^-1 , 



(2.25) 



where we have used the matrix elements of the generators of SU{2), i.e. {tj)mjmj = rrij 5°''^ . 
More generally the exponent on the r.h.s. of Eq.( 2.25 ) can be written as 



dx^m,{x)[A^{x)f^ =2 dx^mj{x)tT[t]A^{x)]. (2.26) 



This gives the path integral representation for Wilson loops defined for SU(2) in the following 
form 



W,{C) = . . ^ / l[Dn,{x) Y^^2j + l)e'^'3fcdx,m,{x)tT[t^^A^ix)]^ 

(2i + i) J ^^(j 

(2.27) 



The exponent contains the gauge field projected onto the third axis in colour space tr[tj74^(x)]. 
Nevertheless, Eq.( p.27 ) differs form Eq.(23) of Ref.[3] by a summation over all values of the 
colour magnetic quantum number rrij of the given irreducible representation j. 

The repeated application of Eq.(^I^) induces that the integrations over the gauge function 
at Xfc are completely independent of the integrations at x^^i. There is no mechanism which 
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leads to gauge functions smoothly varying with Xk (k = 1, . . . n). In this sense the situation is 
opposite to the quantum mechanical path integral. In Quantum Mechanics the integration over 
all paths is restricted by the kinetic term of the Lagrange function. In the semiclassical limit 
fi ^ due to the kinetic term the fluctuations of all trajectories are shrunk to zero around a 
classical trajectory. However, in the case of the integration over gauge functions for the path 
integral representation of the Wilson loop, there is neither a suppression factor nor a semiclassical 
limit like h ^ 0. The key point of the application of Eq.( p.2p and, therefore, the path integral 
representation for Wilson loops is that all integrations over O(xfc) (k = 1, . . . ,n) are completely 
independent and can differ substantially even if the points, where the gauge functions 0,{xk) and 
r2(xfc-i) are defined, are infinitesimally close to each other. 

For the derivation of Eq.(23) of Ref.[3] Diakonov and Petrov have used at an intermediate step 
a regularization drawing an analogy with an axial-symmetric top with moments of inertia I± and 
/|| . Within this regularization the evolution operator representing Wilson loops has been replaced 
by a path integral over dynamical variables of this axial-symmetric top which correspond to 
gauge degrees of freedom of the non-Abelian gauge field. The regularized expression of the 
evolution operator has been obtained in the limit I±, /y — > 0. The moments of inertia have been 
used as parameters like fi ^ 0. Unfortunately, as we show in Sect. 6 the limit I±,I\\ — *■ has 
been evaluated incorrectly. 



3 The SU{N) extension 



The extension of the path integral representation given in Eg. ( 2. 24 ) to SU{N) is rather straight- 
forward and reduces to the evaluation of the character of the matrix C/^(Cxj.xj._i) in the irre- 
ducible representation r of SU{N). The character can be given by [12] 

X[U^{C.,.,^,)] = tr(e^^^="i'^^'^^[C',^-^-;^'']) = 



T.^^.e''^r-^C^,.,,,;A^]^ (3.1) 



rrir 



where Hi (i = 1, . . . , N — 1) are diagonal dr x dr traceless matrices realizing the representation 
of the Cartan subalgebra, i.e. [Hi, H^i] = 0, of the generators of the SU(N) [12]. The sum runs 
over all the weights rrir = {rriri, ■ ■ ■ ,mrN-i) of the irreducible representation r and is the 
multiplicity of the weight rhr and J^fAr = '^r- The components of the vector 
are defined by 

^,[Cx,.,_,;A'']=g I iptHx)], (3.2) 

where we have introduced the notation uj{x) = t°'uj°'{x) = dz ■ A^{x). The functions ip£ [uj{x)] 
are proportional to the roots of the equation det uj{x) — A =0. 

The path integral representation of Wilson loops defined for the irreducible representation r 
of SU{N) reads 

WriC) = Y[ Dnrix) Y: drJ^^(^x^e^9§C^rix)-0Hx)]_ (3 3) 

»' xeC {rnr{x)} 

Let us consider in more details the path integral representation of Wilson loops defined for the 
fundamental representation 3 of SU{3). The character X3[t^3^(C)] is defined as 

X3[[/3^(C)] =tr(e^^i^i[C;^'']+^^2^2[C;A^]) ^g-2«>2[C;^^]/3 
J^,[C;A^]/2V3^i<l>2[C;A^]/6^^-i<^>i[C;A^]/2V3J'^2[C;A^]/6^ (3_4^ 
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where Hi = t^/y/2> and H2 = t^/\/3 [12]. For the representation 3 of SU{3) the equation 
det[a; — A] = takes the form 



- X-truj^{x) - deta;(x) = 0. 



The roots of Eq.( |3.5D read 
A(^) = 



1 



1 



3 



^ ti Lo'^(x) sin - arccos^ 



A(2) 



1 

1 



3 
1 



A(3) 



H — 1= \/tvio'^(x) sin - arccos^ 
V2 V I 3 



1 



- Y tr 0^2(2;) cos ( - arccos^ 



2det 


\ 1 ^^t-tiif'u^^^))' 
tr Ci;2(x) 


2det 


\ 1 ^^t'^tr(t-a;2(^))- 
tr a;2(x) 


2det 


\ 1 ^,rtr(i'^a;2(x))' 

tr 1^2(^2;) 


2det 


1 ^^t'^tr(t-a;2(x))" 
tr ci;2(^2;) 


2det 


i'^tr(t''Lj2(x))' 
1 + 12 / 2/ \ 



In terms of the roots A'^*) {i = 1, 2, 3) the phases $i^2[C'; are defined as 



$i[C;A^] = -gVd (j) y'tr(j2(x) sin ^^arccos^ 

J c 



X ] cos I — arccos 




(3.5) 



(3.6) 



where tra;2(x) = ^ g^y[A^\{x) dx^dxy. Thus, in the fundamental representation 
integral representation for Wilson loops reads 



(3.7) 
3 the path 



- I n [^^^3(x)x 3] (^^i[C;A^]/2^^r^,[C-A^]/& 



+e" 



(3.8) 



where the phases ^>i^2[C;^ ] are given by Eq.(3.7). 



4 Wilson loop for pure gauge field 

As has been pointed out in Ref.[3] the path integral over gauge degrees of freedom representing 
Wilson loops is not of the Feynman type, therefore, it depends explicitly on how one "under- 
stands" it, i.e. how it is discretized and regularized. We would like to emphasize that the 
regularization procedure applied in Ref.[3] has led to an expression for Wilson loops which sup- 
ports the hypothesis of Maximal Abelian Projection [14]. According to this hypothesis only 
Abelian degrees of freedom of non-Abelian gauge fields are responsible for confinement. This 
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is to full extent a dynamical hypothesis. It is quite obvious that such a dynamical hypothesis 
cannot be derived only by means of a regularization procedure. 

In order to show that the problem touched in this paper is not of marginal interest and 
to check if path integral expressions that look differently superficially could actually compute 
the same number we evaluate below explicitly the path integrals representing Wilson loop for 
a pure SU{2) gauge field. As has been stated in Ref.[3] for Wilson loops C a gauge field 
along a given curve can be always written as a "pure gauge" and the derivation of the path 
integral representation for Wilson loops can be provided for the gauge field taken without loss 
of generality in the "pure gauge" form. We would like to show that for the pure SU{2) gauge 
field the path integral representation for Wilson loops suggested in Ref. [3] fails for a correct 
description of Wilson loops. Since a pure gauge field is equivalent to a zero gauge field Wilson 
loops should be unity. 

Of course, any correct path integral representation for Wilson loops should lead to the same 
result. The evaluation of Wilson loops within the path integral representation Eq.( |2.8D is rather 
trivial and transparent. Indeed, we have not corrupted the starting expression for Wilson loops 
(|2.lD by any artificial regularization. Thereby, the general formula (p.8[) evaluated through the 



discretization given by Eqs.(2^) and ( |2.6| ) is completely identical to the original expression 
( p.ll ). The former gives a unit value for Wilson loops defined for an arbitrary contour C and an 
irreducible representation J of SU{2): Wj{C) = 1. 



Let us focus now on the path integral representation suggested in Ref. [3] 



WjiC) = Jl[ Dn{x)e'^''^9^cdx,tr[t^A^{x)]^ 



where all matrices are taken in the irreducible representation J. Following the discretization 
suggested in Ref. [3] we arrive at the expression 

" r 2iJg fr. dx^tilt^A^ix)] 

Wj(C) = lim T] Dn(xk)e ^^^^^+1^^ " ^ ^'^'^ (4.2) 

k=l-' 

Setting A^{x) = d^U{x)U\x)/ig we get 

A'/lix) = ^d^{n{x)U{x)){n{x)U{x))l (4.3) 
By a gauge transformation Q{x)U{x) $7(x) we reduce Eq.([4.1|) to the form 

Wj{C) = / n Dn{x)e'^'^ §cdxt,tT[t^d^n{x)n^x)]_ (4 4^ 

For simplicity we consider Wilson loops in the fundamental representation of SU{2), Wi/2{C)- 
The result can be generalized to any irreducible representation J. 

For the evaluation of the path integral Eq.( |4.4D it is convenient to use a standard s- 
parameterization of Wilson loops C [2]: x^ x^{s), with s E [0, 1] and 2;^(0) = x^{l) = x^. 

The Wilson loop (4.4) reads in the s-parameterization 



Wi/2{C) = J Yl DQ{s) exp J ds ti 



^ dn{s 



0<s<l 



ds 



(4.5) 



10 



The discretized form of the path integral (^]^) is given by 



W,j{C) = hm / TT DVtk expAsfc+i fctr 



k=l 



3 ~ "fc 



+l,k 



= hm / fr £'Jlfce*^[*^^^+^^^'^ = lim / D^nD9-n-iD9-n-2 ■ ■ ■ D^i 

n^oo J n~*oo J J 

fc=l 

^ gtr[t30„f]|,_ J ^2 Jtr[t3j7„„if]|^ J . . . e^At^^2^\] e'^^lt^'^i^l] ^ (4.6) 



where ^n+i = ^i- 

For the subsequent integration over O^, we suggest to use a formula of Ref. [15] modified for 
our case 



(4.7) 



where the coefficients aj (z) are defined by the expansion [15] 



(4.8) 



In the particular case z = 2J and for the fundamental representation J = 1/2 we have z = 1. 
The trace tr[t^C/] in the exponent of the l.h.s. of Eq.( ^.8D should be evaluated for the fundamental 
representation of SU{2). By virtue of the orthogonality relation for characters [9,10,15] 



DU Xj [AU^] Xf [UB] = xj [AB] , 



(4.9) 



where DU is the Haar measure for the SU{2) group, the coefficients aj{z) for j 7^ can be 
determined by [15] 



aj{z) 



j{j + 1) 



- j DU xAt^u^W'^^^'^y 



(4.10) 



We have used here that Xi[(* ) ] = i(j + l)(2i + l)/3. The coefficient ao(-z) is defined by 



a,{z) = j DU e^'^^y^'^y 



(4.11) 



The coefficients aj{z) obey a completeness condition. For its derivation we notice that tr[t^ [/"'"] = 
— tr[t'^C/] which can be easily seen from the standard parameterization of the matrix U in terms 
of an angle (/? and a unit vector n (see Eq.(3.96) of Ref. [10]) 



jj\ ^ e-^V'^-^/2=cos^-i(n-f)sin^, 

2 ^ ^ 2' 

The expansion of e^^^ ^ we represent as follows 

eMt'un^Y.b,{z)xAt'u\ 



(4.12) 



(4.13) 
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Let us show that bj{z) = aj{z). Using the orthogonahty relation ( [4.S| ) we obtain 



boiz) 



j{j + 



(4.14) 



Then, making the change C/^^ — > [/ we get bj{z) = aj{z) by virtue of the invariance of the Haar 
measure DU'^ = DU . 

Thus, taking the product of the expansions ( [4.8| ) and ( 4.13| ) with hj{z) = aj{z) and integrat- 
ing over U we get 



j j' 



aj{z). 



(4.15) 



j>0 



The l.h.s. of Eq.( [4.15| ) is equal to unity due to the relation tr[t^C/^] = — tr[t^C/] and the normal- 
ization of the Haar measure / DU = 1. Therefore, the completeness condition for the coefficients 
aj{z) reads 



]iz) + Y.-jij + l)a]iz) = l. 



(4.16) 



j>0 



The coefficient ao{z) we evaluate below. The evaluation of coefficients aj{z) for an arbitrary j 
is given in the Appendix. For the evaluation of ao{z) one can use, for example, the standard 
parameterization ( 4.12| ) and the definition of the Haar measure DU (see Eq.(3.97) of Ref.[10]) 



1 2 
DU = -— r dQfi dip sin — , 
47r^ 2 



(4.17) 



where dn^i is the uniform measure on the unit sphere S"^ [9]. As a result for a^i^z) we obtain 

,ziT[t^U 



00(2) 



DUe' 



2Ji{z)/z, 



(4.18) 



where Ji{z) is a Bessel function [16]. In the particular case, z = 1, we get ao(l) = o-q = 2Ji(l) = 
0.88 [16]. 

For the integration over 0^, we suppose, first, that n is an even number. Then, integrating 
over iln-i) ^ra-3v • •) ^1 we obtain 



Wii2{C) = Jiirn^ / • • • / DaaDVLn-2D9.n-A ■ ■ ■ DVL2 
X ^(2j„_i + l 



3n-l 



Jn—l 



2i„-i + 1 




Xn[{t'?^2^l] 



(4.19) 
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where we have denoted aj{l) = aj. 

After the integration over r2„, r2„_2,. . ., we arrive at the expression 



2j + l 



hm U + J2(2j + 1) 

j>0 



n 



2j + l 



j>0 



2j + l 



3\n] 



(4.20) 



where we have used that hm. 



by virtue of the relation ao = 2Ji(l) = 0.88 < 1 given 
by Eq.(|418D. 

The expression ( 4.201) is vahd too if n is an odd number. However, in this case Xj[(^^)"] = 
and we get immediately Wi/2{C) = 0. 

In order to estimate Xj[(^^)"] ^or n an even number and at n ^ oo we suggest to apply the 
following procedure 



x,[(t'n= E 

r(n + i) 



m 



m=-j 
ioo 



2TTi 



m=-j 



„ms 



ds 



r(n+l) 7sh(2j + l)- 



2-Ki 



sh — 
2 



(4.21) 



As n — > oo, we can evaluate the integral over s by using the saddle-point approach and get 

x,[(t'n=^i". 

The Wilson loop is then defined by 



^i/2(C^)= li5L5:(2j + l) 



i>0 



2i + l 



lim y(2j + 1) exp 



j>0 



nin 



2j + l 



(4.22) 



By using the completeness condition for the coefficients aj given by Eq.( [4.16| ) we obtain the 
constraint 



2j + l 



< V3(l-a2). 



J 



1 



j + 12j + l 



< 



1 



j + 1 2 j + 1 



< 1. 



(4.23) 



This proves that the Wilson loop Wi/2{C) vanishes in the limit n ^ oo, Wi/2{C) = 0. 

Thus, the Wilson loop Wi/2{C) for an arbitrary contour C and a pure gauge field represented 
by the path integral derived in Ref.[3] vanishes, instead of being equal to unity, Wi/2{C) = 
1. This shows that the path integral representation suggested in Ref.[3] fails for the correct 
description of Wilson loops. 



5 Wilson loop for Z{2) center vortices 

In this Section we evaluate explicitly the path integral (4.1) for Wilson loops pierced by a Z{2) 
center vortex with spatial azimuthal symmetry. Some problems of Z{2) center vortices with 
spatial azimuthal symmetry have been analysed by Diakonov in his recent publication [17] for 
the gauge group SU{2). In this system the main dynamical variable is the azimuthal component 
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of the non-Abehan gauge field A'^{p) (a = 1,2,3) depending only on p, the radius in the 
transversal plane. For a circular Wilson loop in the irreducible representation J one gets 

Wr(p) = ^ V e^2^"i^(/>) = 1 sin[(2J + l)7r^(p)] 

WAP) 2J + 1^^_/ 2J+1 sinMp)] ' ^ ^ 



where p{p) = p^A'^{p)A'^{p). The gauge coupling constant g is included in the definition of the 
gauge field. For Wilson loops in the fundamental representation J = 1/2 we have 

Wi/2(p) = cos[^/i(p)]. (5.2) 

In the case of Z{2) center vortices with spatial azimuthal symmetry and for the fundamental 
representation of SU (2) Eq. (|4.2|) takes the form 



Wi/oip) = lim 17 / DQ^^tT[tHi'^T^p/n)nk+iA^{p)nl + t^nk+inl] 

= lim / . . . / DVLnD9.n-iD9.n-2 ■ ■ ■ D^i 



X e 
X e 
X e 
X e 



tv[t^{i2TTp/n)nnA4p)nl^^ + t^nn^i-i] 

tT[t^i27:p/n)nn-iA^{p)^i^2 + 
tT[t^{i2TTp/n)n2A^{p)n\ + 1^92^1] 

tT[t^{i27:p/n)niA^{p)nl + t^Qinl] ^ (5 3) 



where we have used Cxi^^-i^x^ = 27rp/n, 9{xk) = and i^n+i = f^i- 

For the subsequent evaluation it is convenient to introduce the matrix 

QiA^) = [l + i'^pA4p)). (5.4) 

In terms of Q{A^) the path integral ( |5.3D reads 

Wy2{p) = Jim^l . . . 1 DnnDnn-lDnn-2 ■ ■ ■ D^i etr[^'^«Q(^*)^n-l] 

^ gtr[t30„_iQ(^^)J7],_2] ^ _ _ ^tiit^n^QiA^M] g tr [t3j7iQ(A^)Ot^] ^ ^^^^^ 
The integration over Qk we carry out with the help of Eq.( |4.7| ) taken in the from 

J DVLk e^At'^^k+iQ{A^)^l + Q{A^)^^_-^V^nk] ^ 

= E xA{t'?^k+iQ\A^)9l_,] (5.6) 



and the orthogonality relation ( [4.9[ ). The coefficients aj obey the completeness condition ( [4.16 ) 
with the constraint ( 4.23|) . 

The number n may be both even or odd. Let n be an even number, then integrating over 
^n-i, • •, ^^1 by using Eq.(^.6D we obtain 



Wij2{p) = Jirn^ / • • • / DaaDVLn-2DVLn-A . ■ ■ D9.2 
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X ^(2j„_i + l 

X ^(2j„_3 + l 

in -3 

X ^(2i4 + 1 

J4 



2j„-l + 1 




X 5](2j2 + l) 
The integration over r2„, r2„_2, . . ., gives 
Tyj(p) = lim Y.{2j + 1 



2i„_3 + 1 

2 
2 



n~4i 



(5.7) 



lim (a[J + 5I 

j>0 



lim 7 

r) — ^rm ' ^ 



2j + l 

-1 n 



DnnXj[{tY^nQ''iA^)^i] 

X,[(i')"]x,[Q"(^*)]: 



2j + l 



2j + 1^ 

Xi[(i')"]xi[Q"(^0)], 



(5.8) 



where we have set lim„,_^oo = 0. 

The integration over Qn we have carried out by means of Eq. ( |2.3D . One can easily show that 
Eq.(|5.8|) is also valid for odd n, as well as Eq.( 4.12D . In this case due to the relation Xi[(^^)"] = 
we obtain again Wi/2{p) = 0. 

For even n we should use the relation Xj[(^^)"'] — j" at n — > oo which follows from Eq.( [4.13 ). 
This reduces the r.h.s. of Eq.( |5.7D to the form 



j>Q 



2j + l 



(5.9) 



The evaluation of the character Xj[Q'^iA(f,)] for n — > oo runs as follows 



xAQ'^iA^)] = X.[(l + ^f M0(p))"]-X.[e^2vrM0(p) 



sin[(2i + l)TTfi{p)] 
sin[vr/i(p)] 

Substituting Eq.( |5.10|) in Eq. (|5.9D we obtain 



(5.10) 



^i/2(p)=.lim E 



i>0 



2j + l 



sin[(2j + l)TTfi{p)] 
sin[vr^(p)] 



(5.11) 

Due to the constraint Eq.( 4.23| ) the Wilson loop vanishes in the limit n — > oo, Wi/2{p) = 0. 
Thus, we have shown that the the path integral for Wilson loops suggested in Ref.[3] gives zero 
for a field configuration with a Z{2) center vortex, Wi/2(p) = 0) instead of the correct result 
VFi/2(p) =cos7r/i(p), Eq.(ll). 

We hope that the examples considered in Sect. 4 and 5 demonstrate that the path integral 
representation for Wilson loops derived in Ref.[3] is erroneous. Nevertheless, in Sect. 6 we 
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evaluate explicitly the regularized evolution operator Zjicg{R2, Ri) suggested by Diakonov and 
Petrov for the representation of the Wilson loop in Ref.[3]. We show that this regularized 
evolution operator Z^eg{R2, Ri) has been evaluated incorrectly in Ref.[3]. The correct evaluation 
gives Zjieg{R2, Ri) = which agrees with our results obtained above. 



6 Path integral for the evolution operator Z(i?2,i?i) 

As has been suggested in Ref.[3] the functional Z{R2, Ri) defined by (see Eq.(8) of Ref.[ 



Z(i?2, Ri 



Ri 



R2 / t2 \ 

DR{t) exp liT J Ti (iRRrs)], 



(6.1) 



where R = dR/dt and T = 1/2,1,3/2,... is the colour isospin quantum number, should be 
regularized by the analogy to an axial-symmetric top. The regularized expression has been 
defined in Eq.(9) of Ref.[3] by 



R2 



t2 



ZReg{R2, Rl 



DR{t) exp i 



(6.2) 



Rl 



where Via = iTi[RRTa) are angular velocities of the top, Tq are Pauli matrices a = 1,2,3, I± 
and /|| are the moments of inertia of the top which should be taken to zero. According to the 
prescription of Ref.[3] one should take first the limit /y and then I±_ — > 0. 
For the confirmation of the result, given in Eq.(13) of Ref.[3], 



^Reg(i?2, Rl) = (2T + 1) D^t{R2R\), 



(6.3) 



where D'^{U) is a Wigner rotational matrix in the representation T, the authors of Ref.[3] 
suggested to evaluate the evolution operator ( |6.2| ) explicitly by means of the discretization of 
the path integral over R. The discretized form of the evolution operator Eq. (|6.2| ) is given by 
Eq.(14) of Ref.[3] and readsQ 



N 

ZReg{RN+uRo) = lim M r\ dR 
5 



00 




X exp 



N 

E 

Ln=0 



26 



(TrKTi)2 + (TVKT2)' 



T(TrKT3) 



(6.4) 



where Rn = R{sn) with Sn = ti+n6 and = ^Tr (RnRn^iTa) /S is the discretized analogy of the 
angular velocities [3] and Vn = RnRn^i are the relative orientations of the top at neighbouring 
points. The normalization factor M is determined by 



N+l 



27ri6 V 27ri5 



(6.5) 



(see Eq.(19) of Ref.[3]). According to the prescription of Ref.[3] one should take the limits 5 — > 
and III, /_L — > but keeping the ratios Ii/6, where (i = ||, -L), much greater than unity, Ii/6 » 1. 



^We are using the notations of Ref. [3] 
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Let us rewrite the exponent of the integrand of Eq.( |6.4D in equivalent form 

^+1 , N 



Zueg{RN+i,Ro) = lim 

N ^ oo 

6^0 



X exp 



N 



I± 



2-1116 V 27ri6 



h - 1± 



n dRn 



n=l 



(TrVnT^f- T(TrVnT3) 



n=0 



(6.6) 



26 ' ' 26 

Now let us show that if Vn is a rotation in the fundamental representation of SU{2), so 

{TlVnTaf = -4+{TrVnf. (6.7) 

For this aim, first, recall that 

TT{VnTa) = -TT{V^Ta). (6.8) 

Since is a rotation matrix in the fundamental representation of SU (2) , it can be taken in the 
general standard parameterization given by Eq.( 4.12| ). By virtue of the relation ( |6.8| ) we can 
rewrite (TV VnTa)^ as follows 

(TV VnTaf = -Tr (Kra)TV (V^Ta) = -2 TV ((k - ^ TV Vn) {V^ -^TTVn)) = 
= -2Ti{RnRl^,Rn+iRi) + (TrK)' = -2Trl + (TrK)' = -4 + {TvVnf. (6.9) 
By using the relation Eq.( |6.7D we can recast the r.h.s. of Eq.( |6.6| ) into the form 



ZYieg{RN+i,Ro) = lim 

N ^ oo 

6^0 



2TTi6 V 2TTi6 



N+l 



exp iN{N + 1) 



TV 



Y[dRn 



exp 



n=l 



N 



n=0 \ 



26 



(6.10) 



Now let us proceed to the evaluation of the integrals over i?„ (n = 1, 2, N). For this aim it is 
convenient to rewrite the r.h.s. of Eq. (|6.1C| ) in the following form 



Z^cg{RN+i,Ro) = lim 

iV ^ oo 

6^0 



N+l 



2TTi6 V 2TTi6 



exp iN{N + 1) y 



dR]\f dRj\f^i . . . (ii?2 dRi 



X exp 



26 



(Triijv^Ui)' + {TvRn-iR'j^Y + ... + {TTR2Riy + {Tr R,R, 



?t \2 



t^2 



Ill -/ 



26 



{TrRNR^j^^^Tsf + {Tt RN-iR^^T^f + ... + (Tri^s^-^a)' + (Tr i^ii^Jrg^^ 



-T 



Tr (RnR^n+i^s) + TV (Rn-iRIts) + . . . + TV iR2R\r3) + Tr (i^iiJ^rg) 



(6.11) 
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In the fundamental representation and the parameterization [18] (see Appendix) we have 
TrK = Tr(i?„4^i) = 

(3n Pn+1 ( Oin + ln "n+l + In+l \ 
= I COS COS COS 

2 2 2 2 y 

, „ . /3n . /3n+l ( OLn-ln - 7„+i \ 

+2 sm — sm cos = 

2 2 2 2 y 

( f3n- Pn+1 \ fan- On+l \ fin- In+l \ 

= 2cos ^ j cos^ j cos^ j 

-2 cos ( ^Il±£ll±l\ sin ( ""-""+^'1 sin ( lll^Jll±l] , 
[ 2 ) [ 2 ) [ 2 

Tr (Krs) = TV (RnRl+^r^) = 

„ . Pn f^n+l . f Cln+ln + ln+1 \ 

= —2 1 COS — COS sm 

2 2 2 2 y 

, ^ . . f3n . (3n+l . fan- In On+l - 7n+l \ 

+2 1 sm — sm sm = 

2 2 \^ 2 2 y 

= -2ZCOS ( ^"-^"+^'1 cos ( """""+0 sin f ^""^"+^'| 
V 2 y 1, 2 y 1, 2 y 

„ . fPn + I3n+1 \ . fan- «„+! \ / 7n " 7n+l \ , „^ 

-2^cosl ^ Ism I I cos I I. (6.12) 

The Haar measure Rn is defined by (see Eq.(A.2)): 

DRn = ^ sin Pn dPn dan d'Jn- (6.13) 

Due to the assumption Ii/5 ^ 1, where (i = ||, -L), the integrals over Rn are concentrated around 
unit elements. Expanding Tr (Vn) and Tr (K1T3) around unit elements we get 

TV K = Tr (RnRi+i) =2-^{pn- Pn+lf - ^ (an - an+1 + 7n - 7n+l)^ 
Tr (Vnn) = Tr (RnRl^in) = - i{an- On+l + 7n - 7n+l)- (6.14) 

For the subsequent integration it is convenient to make a change of variables 

an + In 

2 ^ 

an - In an- (6.15) 

The Jacobian of this transformation is equal to unity. After this change of variables ( |6.14D reads 

TV K = Tr (i?„iil+i) = 2 - i (/5„ - Pn+if - {in - ln+l)\ 

TV(Kr3) = TT{RnRi^^n) = -2i{-in-ln+i). (6.16) 

Since both TVy„ and TV(y„r3) do not depend on a„, we can integrate out a„. This changes 
only the Haar measure as follows 

DRn = ^Pndf3nd-in- (6.17) 

47r 
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The integration over /3„ and 7„ we will carry out in the limits — oo < /?„ < oo and —oo < 7n < oo. 
Substituting expansions (6.16) in the integrand of Eq.( |6.1l| ) we obtain 



Z-RcgiRN+i-, Rq 



oo oo 



lim f — —t\ — ^: 

jV — > oo V 27ri() 

5^0 



oo oo 



oo oo 



-oo — oo 



— oo — oo 



~oo — oo 



X exp \i ^ [(/5jv+i - + - f3N~if + . . . + (/^a - + - /3o)'] 

+ ^77? [(7Af+i - iNf + (7Af - 7Af-i)^ + • ■ • + (72 - 71)^ + (71 - 7o)^] 
-2 i T [(7iv+i - 7Ar) + (7Ar - 7Ar-i) + • • • + (72 - 7l) + (7l - 7o)]^ = 

47r / 



p-2iT(7Ar+i -70) 

^ _j. QQ \ 27ri(^ V 27rz(5 
00 00 00 00 



00 00 



00 00 



X / d-iN / d(3N I3n / d7A^-i / d^N-lPN-l ••• / ^72 / d/32 /?2 / dji / 



-00 —00 



—00 —00 



-00 —00 



I± 



X exp ( i ^ [{Pn+1 - Pn? + (Pn - Pn-i? + ... + {p2- Pif + {Pi - Pof] 



+ « ^ KlN+l - iNf + {iN - iN-lf + • • • + (72 - 71)^ + (71 - 7o)^] 



(6.18) 



The integration over 7^ gives 



00 00 00 00 

d^yN j djN-i . . . y ^72 J d'ji 
00 —00 —00 —00 



X exp U ^ [(7Af+i - 7Af)^ + ilN - 7N-i)'^ + • • • + (72 - 71)^ + (71 - 70)^] 



I27ri6 1 l2TTi6 2 l2TTi6N-l l2-Ki5 N I\ 

exp I 



i|i 2y /|| 3 

N 



l2TTi6 



1 



exp I ■ 



N y i|| iv + i 



2{N + 1)5 



ilN+l - 70)' 



/|| J \l N + l 2{N + l)S 

By taking into account the normalization factor the result of the integration over 7^ reads 



(6.19) 



N+l 00 



00 00 



2TTi5 



djN J djN~i ■■■J d-f2 J dji 

00 —00 —00 —00 



X exp ( i — [{-yN+i - Jn^ + (7a^ - iN-i)^ + • • • + (72 - 7i)^ + (71 - 7o)^] 
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" '^^P ( ^ n^Ar ! ilN+1 - 70)^ 



2Tri{N + 1)6 \ 2{N + 1)6 



where we have replaced {N + 1) 6 = t2 — ti = At. The obtained result is exact. By replacing 
/|| — > M, 7Ar+i — > Xfc, 70 — > Xa and At ^ (tb — to) we arrive at the expression for the Green 
function, the evolution operator, of a free particle with a mass M given by Eq.(2.51) of Ref.[6]. 

Thus, after the integration over 7„ the evolution operator ZReg(i?7v+i) Ro) can be written in 
the form 

ZKcg{RN+l-,Ro) = 

= """P (^^(^^+1 -^o)') e-^^^^^^+i -^o)i^[/±,/37V+i,/?o], (6.21) 

where F[I_\_, j3N+i-, Po] is a functional defined by the integrals over 

F[I^,I3n+i,I3o\ = 

^) (^) / / dPN-iPN~i ... J dhP2 J dfiiPi 

^ ^ ^ — oo — oo — oo — oo 

X exp (^i ^ [{(3n+i - Pn? + {(3n - Pn-i? + . . . + {^2 - Pi? + (A - /3o)'] j • (6.22) 

Formally we do not need to evaluate the functional F[I±, [3n+i^ Pq] explicitly. In fact, the 
functional F[I^, (Sn+i, Po] should be a regular function of variables Pn+i and /?o restricted in 
the limit I± — > 0. Therefore, taking the limit /y for the evolution operator ZReg(-RAf+i5 -Ro) 
defined by Eq.( |6.21| ) we get 

ZiR2,Ri) = ^ lim ZReg(i?7v+i,i?o) = 0. (6.23) 

i II ,1 ± — ^0 

This agrees with our results obtained in Sects. 4 and 5. 

Nevertheless, in spite of this very definite result let us proceed to the explicit evaluation 
of the functional F[I±, Pn+i, (3o]. It is convenient to rewrite the integrand of Eq.( |6.22| ) in the 
equivalent form 



F[I±,Pn+i,I3o] = lim 

N ^ oo 

6^0 
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(6.24) 
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After k integrations we get 
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By performing integrations we obtain 



(6.25) 
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(6.26) 



where we have replaced {N + 1) 5 = t2 — ti = /S.t. 

Now we can evaluate the derivatives with respect to ji , j2 , . . . , jn- i,jN- Due to the constraint 
I±/S ^ 1 we can keep only the leading order contributions in powers of I±/6 3> 1. The result 
reads 
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Substituting Eq.( |6.27|) in Eq. (|6.24| ) we obtain the functional /Jtv+i, /3o]: 



F[I^,Pn+i,Po] 
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(6.28) 



In order to understand the behaviour of the functional F[I±, P^^i, Pq] in the hmit ^ oo we 
suggest to evaluate the product 



^ / h AT -I- 1 — 
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(6.29) 



at ^ 1 by using the ("-regularization. In the ^-regularization the evaluation of n[/3jv+i5 /5o] 
runs the following way 
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Thus, the function I1[Pn+i, Po] is defined by 



(6.30) 
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(6.31) 



where e = 2.71828. . .. Due to the constraint I±/5 » 1 the Euler angles Pn+i and Po are less 
than unity and the ratio in Eq.(|6.3l|) is always positive 



Pn+1 in — ^ Poin-^ 
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The functional F[I±, Pn^i, (3q] is then defined by 
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(6.33) 



Thus F[I±, Pjst^i, Po] vanishes in the limit oo. This result retains itself even if we change 

the normalization factor of the evolution operator 
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The renormalized functional /Jat+i, /3o], defined by 
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(6.35) 



vanishes in the limit N ^ oo since the Euler angles Pn+i and Pq are small compared with unity 
due to the constraint I±/6 ^ 1 [3]. The vanishing of the functional F[I±, Pi\i^i, Pq] in the limit 

— > oo agrees with our results obtained in Sects. 4 and 5. 

Substituting Eq.(|05|) in Eq.(|^) we obtain 



ZKcg{RN+l, Ro) — 0. 



(6.36) 



This leads to the vanishing of the evolution operator Z(i?2,i?i) given by Eq.( |6.1| ) or Eq.(8) of 
Ref.[3], Z{R2,Ri) = 0. 

Thus, the evolution operator Z{R2,Ri) suggested in Ref.[3] for the description of Wilson 
loops in terms of path integrals over gauge degrees of freedom is equal to zero identically. This 
agrees with our results obtained in Sects. 4 and 5. As we have shown the vanishing of Z{R2, Ri) 
does not depend on the specific regularization and discretization of the path integral. In fact, 
this is an intrinsic property of the path integral given by Eq. (|6.1| ) that becomes obvious if the 
evaluation is carried out correctly. 



7 Evolution operator Z{R2,Ri) and shift of energy levels of an 
axial— symmetric top 

In this Section we criticize the analysis of the evolution operator Z{R2,Ri) carried out by 
Diakonov and Petrov via the canonical quantization of the axial-symmetric top (see Eq.(12) of 
Ref.[3]). Below we use the notations of Ref.[3]. 
The parallel transport operator 



Waf3{t2,ti) 



x{t2) 



Pexp 



P exp 



A{t) dt 



(7.1) 



J al3 
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where A(t) = A'^{x) dx^/dt is a tangent component of the Yang-Mills field and (a = 1, 2, 3) 
are the generators of SU{2) group in the representation T, has been reduced to the form 

Wo.p{t2M) = Dlp{U{t2)U\t^)), (7.2) 

(see Eq.(5) of Ref. [3]) due to the statement [3]: The potential A{t) along a given curve can he 
always written as a "pure gauge" 

Aa^it) =iDl^{U{t)) ^^D^p{u\t)) (7.3) 

(see Eq.(4) of Ref. [3]). 

By using the parallel transport operator Eq.([i'.2|) the Wilson loop Wt{C) in the representa- 
tion T has been defined by 

Wt{C) = Y.Waa{t2,tl) = Y.Dl^{U{t2)UHti)). (7.4) 



(see Eq.(25) of Ref.[3]). In terms of the evolution operator Z{R2,Ri) given by Eq. (|6.1| ) (see 
Eq.(8) of Ref.[3]) the parallel transport operator Wa/3(t2,ti) has been recast into the form 

W°/(t2,ti) = // dRidR2 (2^' + l)Dll^{U{t2)Rl)D^0{RiUHti))Z{R2,Ri), (7.5) 

'''' T',m 

where the index DP means that the parallel transport operator is taken in the Diakonov-Petrov 
(DP) representation. The Wilson loop Wj^^{C) in the DP-representation reads 

W?^iC) = [[dRidR2 i2T' + l)DlUuit2)Rl)D^MuHti))ZiR2,Ri). (7.6) 

T',m,a 

Of course, if the DP-representation were correct we should get Wj^^{C) = Wt{C), where Wt{C) 
is determined by Eq. (|7.4| ). 

The regularized evolution operator Z^Qg{R2, Ri) given by Eq.( |6.3| ) (see also Eq.(9) of Ref. [3]) 
can be represented in the form of a sum over possible intermediate states, i.e. eigenfunctions of 
the axial-symmetric top 

^Reg(i?2, Ri) = E (2«^ + 1) Di,{R2)DUR\) e-*(*2 - h) Ej^^ (7 7) 

J,m,k 

(see Eq.(12) of Ref. [3]), where Ej^ are the eigenvalues of the Hamiltonian of the axial-symmetric 
top 

J(J + l)-m^ (m-r)2 
Ejm = \ (7.») 

(see Eq.(ll) of Ref. [3]). 

As has been stated in Ref.[3]: // we now take to zero — > (first /y, then I±) we see 
that in the sum (12) only the lowest energy intermediate state survives with m = J = T . The 
resulting phase factor from the lowest energy state can he absorbed in the normalization factor 
in eq.(9) since that corresponds to a shift in the energy scale. 

The statement concerning the possibility to absorb the fluctuating factor exp[— z(t2— ti) T/2I±\ 
in the normalization of the path integral representing the evolution operator is the main one 



allowing the r.h.s. of Eqs.(7^) and (|7.6| ) to escape from the vanishing in the limit /j_ — > 0. 
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In reality such a removal of the fluctuating factor is prohibited since this leads to the change 
of the starting symmetry of the system from SU{2) to U{2). In order to make this more 
transparent we suggest to insert ZReg(-R2) -^i) of Eq.(^) into Eq. (|7.6|) and to express the Wilson 
loop Wj^^{C) in terms of a sum over possible intermediate states, the eigenfunctions of the 
axial-symmetric top. The main idea of this substitution is the following: as the Wilson loop is a 
physical quantity which can be measured, all irrelevant normalization factors should be canceled 
for the evaluation of it. Therefore, if the oscillating factor exp[— i {t2 — ti) T/2I±] can be really 
removed by a renormalization of something, the Wilson loop should not depend on this factor. 



Substituting Eq.( |7.7| ) in Eq.(7.6) and integrating over Ri and R2 we obtain the following 
expansion for the parallel transport operator in the DP~representation 



T' 



w. 



E E DZmt2)UHh))e-'('^-'^)ET'rn_ 



Setting a 



rDPu^ ^ ^ ~ E E 

T' m=-T' 

(3 and summing over a we get the DP-representation for Wilson loops 

T' 



(7.9) 



T'm 



(7.10) 



T' m=-T' 



Due to the definition ([tJ) the r.h.s. of Eq.( fno|) can be rewritten in the form 



T' 



a T' m=-T' 



i{t2 - ti)ET', 



(7.11) 



where Wt'{C) is the Wilson loop in the T' representation determined by Eq.([7.4|). The relation 
( |7.11| ) agrees to some extent with our expansion given by Eq.( |7.lD . 

Following Ref.[3] and taking the limit /y ^ we obtain m = T. This reduces the r.h.s. of 
Eq. (|7.1l| ) to the form 



(C) = E^T'(C)exp 

T' 



i(t2 - tl) 



T'{T' + 1) - T2 

2I1 



Now according to the prescription of Ref. [3] we should take the limit I± 
Ref.[3] and setting T' = T we arrive at the relation 



(7.12) 

0. Following again 



(C) 



Wt{C) exp 



i{t2 - ti) 



T 

2TI 



J2D^^iU{t2)U\t^))exp 



^(t2-t,)^ 



(7.13) 



Thus, Wilson loops in the DP-representation differ from original Wilson loops by the oscillating 
factor exp[i(i2 — ti)T/2Ij_]. The only possibility to remove the oscillating factor exp[— i(i2 — 
ti)T/2I±] is to absorb this phase factors in the matrices U{t2) and U\ti) which describe the 
degrees of freedom of the gauge potential A{t) via the relation (|7.2|) . This yields the changes 



U{t2) 



^7(t2) = C/(t2)e^*2T/2I^^ 
U^ti) = uHti)e-^^^ 



(7.14) 



However, the matrices U{t2) and U^{ti) are now the elements of the U{2) group but not of 
SU{2). Thus, the shift of the energy level of the ground state of the axial-symmetric top 
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suggested by Diakonov and Petrov in order to remove the oscillating factor changes crucially 
the starting symmetry of the theory from SU{2) to U{2). Since the former is not allowed the 
oscillating factor exp[— i(t2 — ii) T/2I^_] cannot be removed. As a result in the limit I_l — > we 
obtain 

W;PP(C) = 0. (7.15) 

The vanishing of Wilson loops in the DP-representation agrees with our results obtained in 
Sects. 4, 5 and 6 and confirms our claim that this path integral representation of Wilson loops 
is incorrect. 



8 The non— Abelian Stokes theorem 

The derivation of the area-law falloff promoted great interests in the non-Abelian Stokes theorem 
expressing the exponent of Wilson loops in terms of a surface integral over the 2-dimensional 
surface S with the boundary C = dS [19] 

irVce'9 fcdx^A^ix) ^^^^^^ , (8.1) 

where Vs is the surface ordering operator [19], da^^ [y) is a 2-dimensional surface element in 
4-dimensional space-time, j; is a current point on the contour C, i.e. x E C, y is a point on the 
surface 5, i.e. y e S, and G^i,{y) = dfj,Ai,{y) - di,A^{y) - ig[Afj_{y), Ay{y)\ is the field strength 
tensor. The procedure for the derivation of the non-Abelian Stokes theorem in the form of 



Eq.(8.1) contains a summation of contributions of closed paths around infinitesimal areas and 
these paths are linked to the reference point x on the contour C via parallel transport operators. 
The existence of closed paths linked to the references point x on the contour C is a necessary 



and a sufficient condition for the derivation of the non- Abelian Stokes theorem Eg. (8.1). 

Due to the absence of closed paths it is rather clear that the path integral representation for 
Wilson loops cannot be applied to the derivation of the non- Abelian Stokes theorem. In fact, 
the evaluation of the path integral over gauge degrees of freedom demands the decomposition 
of the closed contour C into a set of infinitesimal segments which can be never closed. Let us 
prove this statement by assuming the converse. Suppose that by representing the path integral 
over gauge degrees of freedom in the form of the n-dimensional integral ( |2.7| ) we have a closed 
segment. Let the segment Cx^x^.i be closed and the point x' belong to the segment Cxj.Xfc_i, 
x' G Ca;j,Xfc_i- By using Eq.(^^) we can represent the character x\U?'{Cxf.Xk-i)\ by 

= X[^{xi)Ur{Cx^.j,>)Ur{Cxi.j,^_^)VL{xi-i)\ = 

= dr f i?17.(x')x[f^r(a;i)C/r(C,,x')^^J(^')]x[f^r(x')f/r-(a'x,_J^^(Xi-l)] = 



dr / Dnr{x')x[U^{Cx,.')]x[U^{Cx'. (8.2) 



This transforms a (n — l)-dimensional integral with one closed infinitesimal segment into a 
n-dimensional integral without closed segments. Since finally n tends to infinity there is no 
closed segments for the representation for the path integral in the form of a (n — 1) -dimensional 
integral. As this statement is general and valid for any path integral representation of Wilson 
loops, so one can conclude that no further non- Abelian Stokes theorem can be derived within 
any path integral approach to Wilson loops. 
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9 Conclusion 



By using well defined properties of group characters we have shown that the path integral over 
gauge degrees of freedom of the Wilson loop which was used in Eq.(2.13) of Ref.[ll] for a 
lattice evaluation of the average value of Wilson loops can be derived in continuum space-time 
in non-Abelian gauge theories with the gauge group SU{N). The resultant integrand of the 
path integral contains a phase factor which is not projected onto Abelian degrees of freedom of 
non-Abelian gauge fields and differs substantially from the representation given in Ref.[3]. The 
important point of our representation is the summation over all states of the given irreducible 
representation r of SU{N). For example, in SU{2) the phase factor is summed over all values 
of the colourmagnetic quantum number nij of the irreducible representation j of colour charges. 
This contradicts Eq.(23) of Ref.[3], where only term with the highest value of the colourmagnetic 
quantum number nij = j are taken into account and the other 2j terms are lost. This loss is 
caused by an artificial regularization procedure applied in Ref . [3] for the definition of the path 
integral over gauge degrees of freedom. 

As has been stated by Diakonov and Petrov in Ref. [3] the path integral over gauge degrees of 
freedom representing Wilson loops is not of the Feynman type, therefore, it depends explicitly on 
how one "understands" it, i.e. how it is discretized and regularized. In order to understand the 
path integral over gauge degrees of freedom Diakonov and Petrov [3] suggested a regularization 
procedure drawing an analogy between gauge degrees of freedom and dynamical variables of 
the axial-symmetric top with moments of inertia I± and ly. The final expression for the path 
integral of the Wilson loop has been obtained in the limit I±,I\\ — > 0. 

In order to make the incorrectness of this expression more transparent we have evaluated 
the path integral for specific gauge field configurations (i) a pure gauge field and (ii) Z{2) center 
vortices with spatial azimuthal symmetry. The direct evaluation of path integrals representing 
Wilson loops for these gauge field configurations has given the value zero for both cases. These 
results do not agree with the correct values. 

One can show that Eq. (|5Tl| ) can be generalized for any contour of a Wilson loop in SU (2) 

Wy2{C) = JllDn{x)e'9kdx,tr[t'A^ix)]^ 



lim 7 

n— >oo 

j>0 



2i + l 



i2j + l)Wj{C)=0, (9.1) 



where Wj{C) in the r.h.s. is defined by Eq. (|2.1| ) in terms of the path-ordering operator Vc- 
Further, the result (9.1) can be extended to any irreducible representation of SU{2). Thus, we 
argue that the path integral suggested in Ref. [3] to represent Wilson loops is identically zero for 
Wilson loops independent on the gauge field configuration, the shape of the contour C and the 
irreducible representation of SU{2). 

This statement we have supported by a direct evaluation of the evolution operator Z^eg{R2, Ri] 
defined by Eq.(14) of Ref. [3], representing the assumption by Diakonov and Petrov for Wilson 
loops in terms of the path integral over gauge degrees of freedom. As we have shown in Sect. 6 the 
regularized evolution operator ZYieg{R2, Ri), evaluated correctly, is equal to zero. This agrees 
with our results obtained in Sects. 4 and 5. In Sect. 7 we have shown that the removal of the 
oscillating factor from the evolution operator suggested in Ref. [3] via a shift of energy levels of 
the axial-symmetric top is prohibited. Such a shift of energy levels leads to a change of the 
starting symmetry of the system from SU{2) to U{2). By virtue of the oscillating factor the 
Wilson loop vanishes in the limit — > in agreement with our results in Sects. 4, 5 and 6. 

We hope that the considerations in Sects. 4-7 are more than enough to persuade even the 
most distrustful reader that the path integral representation for the Wilson derived by Diakonov 



28 



and Petrov by means of special regularization and understanding of the path integral over gauge 
degrees of freedom is erroneous. 

The use of an erroneous path integral representation for Wilson loops in Ref.[7] has led to 
the conclusion that for large distances the average value of Wilson loops shows area-law falloff 
for any irreducible representation r of SU{N). Unfortunately, this result is not supported by 
numerical simulations of lattice QCD [8]. At large distances, colour charges with non-zero 
A'^-ality have string tensions of the corresponding fundamental representation, whereas colour 
charges with zero A^-ality are screened by gluons and cannot form a string. Therefore, the result 
obtained in Ref . [7] cannot be considered as a new check of confinement in lattice calculations as 
has been argued by the authors of Ref. [7]. 

We would like to accentuate that the problem we have touched in this paper is not of marginal 
interest and a path integral, if derived by means of an unjustified regularization procedure, would 
hardly compute the same physical number as the correct one. We argue that no regularization 
procedure can lead to specific dynamical constraints. In fact, the regularization procedure 
drawing the analogy with the axial-symmetric top has led to the result supporting the hypothesis 
of Maximal Abelian Projection pointed out by 't Hooft [14]. Any proof of this to full extent 
dynamical hypothesis through a regularization procedure and through specific understanding of 
the path integral should have seemed dubious and suspicious. 

Finally, we have shown that within any path integral representation for Wilson loops in 
terms of gauge degrees of freedom no non-Abelian Stokes theorem in addition to Eq.( ^.l[ ) can 
be derived. Indeed, the Stokes theorem replaces a line integral over a closed contour by a surface 
integral with the closed contour as the boundary of a surface. However, approximating the path 
integral by an n-dimensional integral at n ^ oo there are no closed paths linking two adjacent 
points along Wilson loops. Thereby, the line integrals over these open paths cannot be replaced 
by surface integrals. Thus, we argue that any non-Abelian Stokes theorem can be derived only 
within the definition of Wilson loops through the path ordering procedure ( |8.1| ). Of course, one 
can represent the surface— ordering operator Vs in Eq. ( ^.ID in terms of a path integral over gauge 
degrees of freedom, but this should not be a new non-Abelian Stokes theorem in comparison 
with the old one given by Eq. ( |8.1D . That is why the claims of Ref. [3-5] concerning new versions 
of the non-Abelian Stokes theorems derived within path integral representations for Wilson 
loops seem unjustified. 

Discussions with Jan Thomassen are appreciated. 
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Appendix. Coefficients aj{z) 



In this Appendix we evaluate the coefficients aj{z) of the expansion (4.8) and show that the 
completeness condition given by the series ( [4. IS] ) converges to unity. 
It is convenient to rewrite Eq.(|4.8|) as follows 



+ai(z) xiit'^U] + 03/2(2) X3/2[t'^U] + 02(2) X2[t'^U] + 05/2(2) X5/2[t'^U] + ... (A.l) 

For the evaluation of the coefficients 00(2), 0,1/2(2), ai{z), 03/2(2), 02(2) and 05/2(2) we would 
make use the parameterization of the matrices U in terms of the Euler angles a, (5 and 7 ranging 
over the regions > a > 27r, > /? > vr and > 7 > 27r, respectively [18]. Then, the Haar 
measure DU normalized to unity is defined by [18] 



DU = — ^ sin j3 df3 da d'j. 



In the notation of Ref.[18] the characters are given by 



Xi/2[t U] = -icos - 



(3 . a + 7 



sm 



{A.2) 



X3/2[t'C/] 





Xi[t^U] = -2i cos^^ sin (a + 7), 

10 ■ 3 /3 . a + 7 2 « + 7 , o . /3 . a + 7 

— IZi cos — sm cos \- 2i cos — sm , 

2 2 2 2 2' 

X2[t^U] = —8i cos^ — sin (a + 7) cos(a + 7) 



+ i (1 + cos /?) (1 — 2 cos 13) sin (a + 7) , 



X5/2[i'f/] = -80i cos^l [ sin ^ 



a + 7 „ . 3 « + 7 , . 5 " + 7 
- 2 sm h sm 



-i I 15 cos^ ^ - 12 cos^ ^ 



Q ■ « + 7 , . 3 Q + 7 
3 sm 4 sm 



- i f 10 cos^ ^ - 12 cos^ ^ + 3 cos ^ 



sm ■ 



a + 7 



(A.S) 

The coefficients 00(2), 01/2(2), 01(2), 03/2(2), 02(2) and 05/2(2) are defined by the integrals 
( pop and ( pTD 



27r 27r tt 



O0l2j 



(3 . a + 7 

— z 2 cos — sm 

dPsmf3e 2 2 , 



2^ 



000 
01/2(2) =4| I)[/xi/2[t'f/^]e^^V2[*'f^] 

27r 27r tt 

da I dj I dp sin /3 cos — sin ^^-^ ^ 



000 



/? . a + 7 

— z 2 cos — sm 

e 2 2. 
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2n 2n n 



J da J dj J dp sin f3 





cos^ ^ sin {a + 7) 



/3 . a + 7 
— z z cos — sm — - — 
e 2 2, 



27r 27r 7r 



a3/2(^) = 5 / 1)^7 X3/2[i't/^]e^ ^1/2 ^] = _1_ y da | ^7 / ^^Z? sin /3 



000 



0/3 . a + 7 oQ^ + T P ■ " + 7 

cos — sin — - — cos — cos — sm 

2 2 2 2 2 



P . a + 7 
—I z cos — sm — - — 
? 2 2 . 



27r 27r tt 



a2(^) = ^/ Z^C/x2[i't/^]e^^V2[*'^] = ^ J da J dj J dp sin P 



000 



4/3 



8 cos — sin (a + 7) cos(a + 7) — (1 + cos /?) (1 — 2 cos /?) sin (a + 7) 

/? . a + 7 

—I z cos — sm 

xe 2 2 , 

27r 27r tt 

a5/2(^) = ^/ DUxs/2[t'U^]e'^^/'^''^^ = 7^ / / / ^ 

000 

5/3/. a + 7 „.s" + 7 .';a + 7\ 
80 cos*" - sm - 2 sm"* + sm^ ' 



+ (15 cos5?-12 cos^^ 



o ■ " + 7 , . 
, , 3 sm 4 sm' 

2 \ 2 2 



2 J 

3 a + 7' 



+ ( 10 cos^ ^ - 12 cos^ TT + 3 cos ^ I sin 



2 



P \ a + 7 
2 



P . a + 7 
— z z cos — sm — - — 
= 2 2 . 



The integration over angle variables gives 

ao{z) = 2 ^ , ai/2(z) = 8 ^ , ai(z) = , 

16 Mz) / ^ 72 J6(z) 

03/2(2) = y a2(z) = , 05/2(2) = — . . . 

One can show that for an arbitrary j the coefficient aj (z) is defined by 

im^^ImM . , = 1/2.3/2,5/2..... 

I , j = l,2,... 

This implies the following general formula 

/z)[/x,[i^t/t]e^^V2l^'f^T = / (2J + 1)^^^^ , j = 1/2, 3/2, 5/2,.. 
J ' I , j = l,2,... 
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The completeness condition ( [4.16 ) for arbitrary z takes the form 

^ -^(-^ + ^) ^ 

The numerical values of ao(^), 0,1/2 (-2), 03/2 (-z) and 05/2(^2) for the particular case z = \ read 

ao(l)=0.88 , ai/2(l) = 0.92 , 03/2(1) = 0.01, , 05/2(1) = 4 x 10"^ (^.9) 
For the completeness condition ( [4.16D we obtain at z = 1 

°° 1 1 5 35 

«0(1)+ E ;Ti(j + l)o,'(l)=og(l) + -0?/2(l) + -o2/2(l) + -o2/2(l) + ...= 
i=l/2"^ ^ ^ 

= 0.774 + 0.212 + 1.25 x 10"^ + 5 x 10"^ + . . . = 0.986 + . . . (^.10) 
Thus, this series converges slowly to unity. 
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